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Abstract 

Electroweak one-loop calculations for production of top-quark pairs at colliders 
are completed by providing the missing QED type contributions from real and 
virtual photons, where also effects from interference between QED and QCD 
contributions have to be taken into account. Moreover, photon-induced tt pro- 
duction is included as another partonic channel. 



1 Introduction 



Experimental investigations of the top quark at the Tevatron have significantly con- 
tributed to precision tests of the Standard Model (SM) since the top discovery in 
1995 [1, 2]. The top quark mass is an important parameter within the SM and its pre- 
cise knowledge is an essential ingredient to constrain the mass of the Higgs boson [3] . 
Besides the top mass, the measurement of the top-pair production cross section is an 
important test of the SM, and possible observation of deviations from the SM predic- 
tions could indicate new, non-standard, contributions. Moreover, precise knowledge 
of the SM processes as a main source of background is crucial in direct searches for 
potential new physics beyond the SM. 

At the Tevatron, the dominant production mechanism is the annihilation of quark- 
antiquark pairs q + q ^ t + t, wheras at LHC energies, tt production proceeds mainly 
through gluon fusion, g + g ^ t + t. In lowest order, the ti production cross section 
in hadronic collisions is of 0{a1) and was calculated in [4]. The corresponding lowest- 
order electroweak contributions of 0{a^) to the Drell-Yan annihilation process via 
7- and Z-exchange are very small, contributing less than 1% at the partonic level [5], 
and are thus negligible. Accordingly, the main higher order contributions arise from 
QCD. Cross sections and distributions including QCD effects of 0{aD were computed 
in [6, 7], and an inspection of the QCD effects close to the production threshold was 
performed in [8]. Including the resummation of large logarithmic QCD contributions 
in the threshold region improves the perturbative calculation and was done in [9, 10, 
11, 12, 13]. The prediction for the tt production cross section currently used at the 
Tevatron is based on the studies in [14], which include the next-to-leading-order (NLO) 
contributions and the resummation of soft logarithms (NLL). In [15], also the next-to- 
next-to-leading-order (NNLO) soft-gluon corrections were taken into account, extended 
to NNNLO in [16]. 

Prom the electroweak (EW) side, the EW one-loop corrections to the QCD-based 
lowest order calculations, which are of 0{aal), were investigated first in [17] for the 
subclass of the infrared-free non-photonic contributions, i.e. those loop contributions 
without virtual photons. They are of special interest due to the large Yukawa coupling 
of the top quark to the Higgs boson. However, they have little impact within the SM, 
about 1% of the lowest-order cross section for the Tevatron, and not more than 3% for 
the LHC [17, 18]. In these calculations contributions including the interference of QCD 
and EW interactions were neglected. A study of the non-photonic EW corrections with 
the gluon-Z interference effects was done more recently in [19, 20, 21, 22]. 

Still, a subset of the full EW corrections, corresponding to the QED corrections 
with real and virtual photons, was not included in the previous calculations. In this 
paper we close this gap and present the calculation of the missing QED subset, thus 
making the SM prediction at the one-loop level complete. 

It is worth to mention also several studies within specific extensions of the SM, 
comprising calculations of the Yukawa one-loop corrections within the General 2-Higgs- 
Doublet Model (G2HDM) for Tevatron [23] and LHC [24]. Also, the SUSY-QCD 0{a^) 
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contributions were calculated for Tevatron [25, 26] and LHC [27], and for both [28]. 
The SUSY-EW corrections have been examined: partial calculations relevant for the 
Tevatron were done in [26], and a more complete description at 0{aag) within the 
G2HDM and MSSM with numerical results for Tevatron and LHC, was presented 
in [29, 30]. 

In the following, we provide the QED corrections to top pair production and also 
the effects arising from interference of QCD and QED interactions that occurs at one- 
loop order. Moreover, at this order the distribution of photons inside the proton has 
to be taken into account, adding photon- induced top production as another partonic 
channel at NLO. In the end, wc present numerical results for both Tevatron and LHC. 

Although our calculations arc performed in the frame of the SM, they are also valid 
for extensions of the SM, as e.g. the MSSM and G2HDM. 

2 Structure of the QED contributions 

The QED contributions can be treated as a separate subclass at the electroweak one- 
loop level. They consist of virtual and real photon contributions, according to the 
topology of photonic insertions in the lowest-order graphs. Both real and virtual photon 
terms have to be combined in order to obtain a consistent, infrared (IR) finite result. 

2.1 Virtual corrections 

The virtual QED corrections consist of loop contributions with virtual photons. They 
can be described by the matrix elements dM.a, o, = gg, qq, for both partonic subpro- 
cesses separately. Contracting these quantities with the Born matrix elements 
yields the one-loop contribution to the differential cross sections at the partonic level 
of the order of 0{aa'^), after spin and colour summation. 



where t and s are the usual Mandelstam variables. The explicit expressions of Ai'^ 
are given in [17], where the non-photonic electroweak corrections have been studied. 
Throughout this paper we closely follow the notation of [17]. 

The virtual QED corrections of 0{aa'^), contributing to SMa, can be classified 
according to self-energy, vertex, and box corrections, depicted by Feynman diagrams 
in Figs. 1 and 2 for the two partonic processes of qq annihilation and gluon fusion. They 
were treated with the help of the FeynArts [31], FormCalc [32], and LoopTools [33] 
packages, based on techniques from [34, 35], which were further refined for 4-point 
integrals in [36, 37]. The analytical expressions for the matrix elements are close to 
those in [17] and [38]. 

The whole set of QED loop diagrams is gauge invariant and UV-finitc after taking 
into account the counter terms for the gqq- vertex, gti-yeitex and top quark self-energy. 
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Figure 1: Virtual QED 0{aa'^) contributions to gluon fusion (li-channel diagrams are 
not explicitly shown). Crossed lines and vertices denote counter term insertions. 
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Figure 2: Virtual QED 0{aal) contributions to qq annihilation 



UV singularities in the sum of vertex functions and corresponding counter terms with 
quark field and mass renormalization constants cancel, hence, no coupling constant 
renormalization is needed. 

To obtain finite vertices and propagators, it is thus sufficient to perform field and 
mass renormalization for the quarks. In case of top quarks, the substitution 



rrit 



1 + -5Ztj , 
rrit - Smt , 



(2) 



yields the counter term for the g^tt- vertex, 5A^, and for the top self-energy, 5E, as 
follows: 
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i5K^ = -igsT^^n SZt , (3) 

— X — ^ i5E = i{i)8Zt-mt8Zt + 5mt) . (4) 

The field renormalization constant 5Zt as well as the mass counter term 5mt are fixed 
by renormalization conditions, for which we choose the on-shell scheme. They are 
imposed on the renormalized top quark self-energy E = S + 5E, with 

rv 

E(p) = -[^E^(p2)+m,S5(p^)] , (5) 
corresponding to the QED-like part of the unrenormalized top quark self-energy, 



(6) 




in the following way: 

(i) The pole of the top quark propagator is kept at mt and thus defines the on-shell 
mass: 

Re E (i^ = m^) = , yielding ^ = -(Ey + E5) | . (7) 

(ii) The residue of the top quark propagator is unity, yielding the field renormalization 
constant by (real parts only) 

SZt = -T.y{ml) - 2ml ^i^v + • (8) 

The renormalization constants for initial-state quarks (q) are determined analo- 
gously, substituting mt — > m^. To obtain the counter terms for the initial state gqq 
vertices, only the renormalization of the quark fields is necessary. Light quark masses 
are only kept where they are necessary to regularize coUinear divergences, which appear 
as double logarithms ln^{s/m'^) and single logarithms ln(s/m^). 

As a consequence of the null photon mass, the virtual QED corrections are IR 
divergent. The photonic IR singularities can be regularized by introducing a fictitious 
photon mass A. 

A specific peculiarity of the QED corrections are the 0{aas) box contributions 
shown in Fig. 2, which contain besides photons also gluons in the loop. As a conse- 
quence, further IR singularities related to the gluons emerge from the loop integrals. 
Since the gluons appear quite similar to the photons in the box graphs, it is possible 
to perform the regularization by a gluon mass as well. For simplicity, we use the same 
regularization parameter A. 
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Figure 3: Real QED ^(aag) contributions of photon bremsstrahlung to the gg fusion 
(u-channel diagrams are not exphcitly shown). 




Figure 4: Real QED 0{aag) contributions of photon bremsstrahlung to the qq anni- 
hilation. 

2.2 Real corrections 

According to the Bloch-Nordsieck theorem [39], the IR singularities in the virtual 
corrections cancel against their counterparts from the real photon contributions after 
integration over the photon phase space. Therefore we have to include all contributions 
of the real photon radiation off the external particles to obtain an IR finite result. The 
corresponding diagrams are shown in Figs. 3 and 4. 

Moreover, we have to include also gluon bremsstrahlung to compensate the IR 
singularities related to the gluons in the box graphs of Fig. 2. They consist of two 
types of diagrams: gluon radiation off the QED-mediated and off the QCD-mediated 
qq annihilation, as depicted in Fig. 5. 

At 0{aa'^), it is the interference of these two classes of diagrams that is required, 
yielding a new type of QED-QCD interference. Still, not all of the interference terms 
contribute. Owing to the color structure, only the interference of the initial and final 
state gluon radiation graphs is non-zero, yielding the structure required to cancel the 
IR singular parts in the box corrections of Fig. 2. Nevertheless, the cancellation is not 
yet complete. The missing piece is the pure QCD box correction interfering with the 
QED qq annihilation Born level diagram, as displayed in Fig. 6, which gives another 
non-zero contribution of the same order. Only after combining all these various parts 
the 0{aa'^) result is IR finite. 

In Fig. 6, only the photon-mediated Born-level diagram is shown. In principle, also 
the interference of the QCD box and Z-boson exchange tree-level diagram has to be 
taken into account. This contribution belongs to the IR-singular gluon-Z corrections, 
which also contain the gluon-Z box graphs and gluon bremsstrahlung off Z-mediated 
tree-level diagrams. The IR-singular structure of these contributions is simplified by the 
fact that there are no IR-singularities related to the Z-boson. The gluon-Z interference 
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Figure 5: Gluon bremsstrahlung from QED-mediated (upper row) and QCD- mediated 
(second row) Born diagrams, contributing at 0{aa:^) through interference. 




Figure 6: 0{aaf) contribution to qq annihilation via QCD box diagrams (crossed 
diagram not exphcitly shown) interfering with the QED Born diagram. 

effects were neglected in the original study of non-photonic EW corrections performed 
in [17]. They have been investigated recently in [19, 20, 21, 22]. 

For completeness of the NLO QED effects, photon radiation off the off-shell top 
quarks in the gg fusion subprocess (Fig. 7) has to be considered as well. These effects 
arc, however, one order of magnitude smaller than the other terms and hence are less 
important for numerical studies. 

Technically, for the phase space integration of real photon/gluon radiation, we 
apply the phase space slicing method (see appendix) taking advantage of its univer- 
sality in handling both inclusive and non-inclusive quantities. The dipole subtraction 
method [40], originally proposed for QCD [41], was used to verify numerical results 
obtained with the slicing method at the partonic level. 




Figure 7: IR- finite 0{aa^) bremsstrahlung contributions. 
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Figure 8: Feynman diagrams for photon induced tt production at lowest order. 
2.3 Photon- induced tt production 

In addition to the previously mentioned NLO QED contributions we also have to 
inspect the photon-induced production channels. These comprise at lowest order the 
gluon-photon fusion amplitudes illustrated in Fig. 8. 

In general, photon-induced partonic processes vanish at the hadronic level unless 
the NLO QED effects are taken into account. A direct consequence of including these 
effects into the evolution of parton distribution functions (PDFs) is the non-zero photon 
density in the proton, which leads to photon-induced contributions at the hadronic level 
by convoluting the photon-induced partonic cross sections with the PDFs at NLO QED. 
Since the photon distribution function is of order a they are formally not of the same 
overall order as the other NLO QED contributions. Numerically, however, they turn 
out to be sizeable, and we therefore include them in our discussion. 

As the PDFs at NLO QED have become available only recently [42], the photon- 
induced hadronic processes have not yet been investigated. Here we present the first 
study of these effects on the top pair production. 

3 Hadronic cross section for pp, pp — > tiX 

For obtaining the hadronic cross section we have to convolute the various partonic 
cross sections with the corresponding parton densities and sum over all contributing 
channels, adding up contributions of the non-radiative and radiative processes. As 
already mentioned, only the sum of all virtual and real corrections is IR finite. Final 
step is the factorization of the remaining mass singularities. 

3.1 Mass factorization 

The mass-singular logarithmic terms proportional to In niq are not canceled in the sum 
of virtual and real corrections. They originate from coUinear photon emission off the 
incoming light quarks. In analogy to the factorization of coUinear gluon contributions, 

they have to be absorbed into the parton densities. 

This can be formally achieved by replacing the bare quark distributions qi{x) for 
each flavor by the appropriate scale dependent distributions qi{x,Q'^) in the following 
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way, according to [43] (with mj = m^.), 
qi{x,Q'^) = qi{x) 



+ ^Qhi^) I 1 - In 5, - ln2 + (^In 6, + In 

+ / —<l^{-)7^Q^S-. 



z \zJ 27r^^ 1 1-;^ ymlil-zf 

1 + ^2 

+ \Fch ) , (9) 



involving a soft-energy cut 5s and the functions 

27r2 

/,+, = 9 + — + 31n5,-21n2 5, , (10) 



1 + ^2 3 

/c= In - + 2^ + 3. 11 

1-z V / 2(1-^) ^ ' 

The expressions for the PDFs are given in both DIS and MS factorization schemes, 
which corresponds to Ai^c" = 1 and \pc = respectively. 

For a consistent treatment of the collinear singularities at 0(«), it is necessary 
to use an appropriate set of PDFs that was extracted from the data and evolved by 
DGLAP equations with the NLO QED effects included. Otherwise it would lead to an 
overestimation of the scale dependence. We use the PDFs from the MRST collaboration 
[42] which were determined at NLO QCD and NLO QED. The authors do not explicitly 
state which factorization scheme is relevant for NLO QED. We follow the reasoning 
given in [44] and use the DIS scheme in our calculation. For the numerical evaluation 
the factorization scale is set to Q = = 2mt. After performing the factorization of 
mass singularities, the results become free of the quark- mass logarithms. The scale 
dependence cannot be checked in a consistent way owing to the NLO QCD effects in 
the parton densities, as these are not included in our calculation. For this reason we 
do not present a study of the scale dependence. 

3.2 Integrated hadronic cross sections 

The observable hadronic cross section is obtained by convoluting the short distance 
partonic cross sections dgg, dq-q, a^g with the universal parton distribution functions 
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for quarks, gluons, and photons. For colliding hadrons A and B carrying the momenta 
Pa and Pb, with S = {Pa + PbY, the hadronic cross section can be expressed as 



»i 



■■•t 
s 



(12) 



with T — s/S, and s the partonic center-of-mass energy squared. The partonic cross 
sections agg and ^gg{s) include the virtual photon and gluon loop contributions as well 
as the real photon and gluon bremsstrahlung terms, as described in section 2; the 
photon-induced partonic cross section a^g is of lowest order. 
The parton luminosities are defined as follows, 

^ rlr \- T 1 

^in{r) = j-^ - $^/^(x,M$n/B(-,M + (l--2) , (13) 

with the parton distributions inside A, ^m/A — Qi, 9, 7, and the factorization scale, /ip- 
3.3 Differential hadronic cross sections 

In addition to the integrated hadronic cross section it is convenient to define hadronic 
cross sections differential in one or more parameters. Typically, the variables are chosen 
to be Lorentz invariant quantities or quantities with simple transformation properties. 
In our study we consider the differential cross sections with respect to the invariant 
mass of the ti pair (inclusive) and with respect to the transverse momentum of top 
quark. 

The invariant mass distribution of the hadronic cross section has the following form, 

|m,n| 

where the sum extends over the various partons m,n in the initial state. 

The differential hadronic cross section with respect to the transverse momentum of 
the top quark, p^, can be written for the partonic 2^2 processes as follows, 

dcr , „, ^ di da ,^ 

drC{T)- r(t,s) (15) 



dpT J TO dt 

(dropping parton indices and summation). The lower limit for the r-integration is 
thereby dependent on px, 

f„ = i(!r!?+ii) . (16, 

Real photonic/gluonic corrections belong to three-particle final states. Therefore, 

we also need the differential hadronic cross section with respect to the transverse mo- 
mentum of the top quark, in each of the 2 — > 3 parton processes. Expressed in terms 
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of variables of the parton CMS, it can be written in the following way, 



where ki and are the 4- momenta of the top quark and the photon(gluon). 6 is the 
angle between ki and the beam axis, given by Pa, and 03 is the azimuthal angle of ^3 
with respect to ki as polar axis. The threshold for the r integration corresponds to 

foS = ( x/mi +pi + \/ [mt + X)^ + pi ] , (1^ 



with the IR mass regulator A. 

4 Numerical results 

In the following we present numerical results for the total hadronic cross section as 
well as for the distributions with respect to the invariant mass of the tt pair and the 
transverse momentum of the top quark. 

For the identification of ti pairs and event reconstruction it is necessary to apply 
kinematical cuts, such as cuts to the transverse momentum and the pseudorapidity 77 
of the t and i. In the case of the LHC, the cuts applied are as follows, 

Pt > 100 GeV and I77I < 2.5 . (19) 

For the Tevatron, the cuts are chosen according to 

PT>25GeV and |r/|<2.5. (20) 

In Table 6.2 and 6.3, we present the numerical results for the integrated hadronic 
cross sections at the LHC and at the Tevatron, respectively. The values of a are 
listed for the Born level and for the NLO QED corrections. The contributions of all 
production channels are shown separately, as well as combined to the total correction. 

At the LHC, the largest correction comes from the photon-gluon production channel 
at NLO. It has the same sign as the contribution to the gg fusion which is the dominant 
ti production channel at the LHC. However, the contributions to qq annihilation have 
opposite signs which leads to a reduction of the overall NLO QED correction. In total, 
the relative correction is about 1% and is slightly reduced if the cuts are applied. 

At the Tevatron, the largest contributions to the total hadronic cross section come 
from the uu subprocess. The photon-gluon subprocess yields the second largest con- 
tribution, but with opposite sign. In total, the relative correction can amount to 2.5%, 
including cuts. 

For illustration of the numerical impact of the NLO QED corrections on the distri- 
butions, we introduce the relative correction S, defined as 
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Table 1: Integrated hadronic cross section for tt production at the LHC, at NLO QED 
in different production subprocesses, without and with cuts. 



Process Utot without cuts [pb] Utot with cuts [pb] 





Born 


correction 


Born 


correction 


uu 


34.25 


-1.41 


18.64 


-0.770 


dd 


21.61 


-0.228 


11.54 


-1.68 


ss 


4.682 


-0.0410 


2.253 


-0.0304 


cc 


2.075 


-0.0762 


0.9630 


-0.0446 


99 


407.8 


2.08 


213.6 


0.524 


91 




4.45 




2.29 


VP 


470.4 


4.78 


247.0 


1.80 



Table 2: Integrated hadronic cross section for tt production at the Tevatron, at NLO 
QED in different production subprocesses, without and with cuts. 



Process 


Ctot without cuts [pb] 




Ctot with cuts [pb] 




Born 


correction 




Born 


correction 




uu 


3.411 


-0.117 




3.189 


-0.118 




dd 


0.5855 


-2.89x10" 


-3 


0.5432 


-2.91x10" 


-3 


ss 


8.063x10-3 


-1.21x10" 


-5 


7.343x10-3 


-1.79x10" 


-5 


cc 


2.044x10-3 


-5.06x10- 


-5 


1.857x10-3 


-5.00x10- 


-5 


99 


0.4128 


3.17x10- 


-3 


0.3803 


2.69x10- 


-3 


91 




0.0154 






0.0143 




PP 


1.120 


-0.102 




1.121 


-0.101 





with cross section at NLO, (io-NLO; and the Born cross section daB- 

In Fig. 9 the and distributions are shown (left), as well as the relative QED 
corrections (right), for the gg and qq parton channel at the LHC. The effect of the 
NLO QED corrections in the dominant gg fusion channel is rather small, less than 1% 
over most of the px range and also over most of the \/§ range. Differently from the gg 
channel, the NLO contributions for qq annihilation arc negative over the whole px and 
range, reaching the 5% level for > 400 GeV and > 1200 GeV. They further 
grow in size with increasing px and \/§ and for very high px the qq channel starts to 
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dominate over the gg fusion. 

In case of the Tevatron [Fig. 10], the qq annihilation dominates over the gg fusion 
(left). The impact of 0{aa'^) corrections on both channels is similar to the LHC. 
Again, in the gg fusion, the relative correction S is smaller than 1% for most of the 
Pt and ^/I ranges, except for the low p-j and threshold regions where it reaches about 
2% (up right). In the qq annihilation channel, the relative corrections are negative 
and at a few per cent level already near the threshold. They grow further in size with 
mcreasmg px and Vs. The 5% level is reached for px > 350 GeV and V§ > 900 GeV. 

As previously discussed, also the photon-induced processes represent contributions 
of NLO in QED, owing to higher order effects included in the PDFs. In Figs. 11 and 12, 
we show the photon-gluon contribution to tt production, in comparison with the NLO 
terms in the gg and qq channels. For the LHC [Fig. 11], the photon- induced contribu- 
tion is larger than the corrections to both Born processes, a consequence of the fact that 
the combination of gluon and photon parton densities can become substantially large. 
Since (77 hadronic cross section is of the same sign as the NLO contributions to the gg 
fusion channel, its presence enhances the size of the overall NLO QED contributions. 

The situation is different at the Tevatron [Fig. 12], where the qq annihilation channel 
dominates. Still, the photon-induced contribution is larger in size than the NLO QED 
corrections to the gg fusion channel. However, as a consequence of opposite signs, these 
two tend to reduce the contribution from the qq channel. 

Finally, we show the combination of the partial results for all production subpro- 
cesses, including the photon-gluon channel, in Fig. 13 for the LHC and in Fig. 14 for 
the Tevatron. As a consequence of the dominant gg production channel at the LHC, 
the total NLO QED corrections in the distribution are positive and at the level of 
about 1%. After applying the cuts and in the px distribution the corrections become 
negative and tend to increase in size with px and v^I. This is caused by the logarithmic 
final state radiation contribution which is not canceled by the virtual corrections in 
case of non-inclusive quantities. 

At the Tevatron, with the dominating qq channels, the total NLO QED corrections 
have a larger impact owing to the subtleties of the QED-QCD interference, which is not 
present at the Born-level. They are negative and in size of several per cent, becoming 
larger with increasing px and \/§, and are slightly enhanced by the application of cuts. 

5 Conclusion 

We have provided the last missing item for a complete EW one-loop calculation for 
ti production at hadron colliders. The NLO QED contributions form, together with 
the non-photonic EW contributions, the complete EW corrections to ti production at 
the one-loop level. For consistency and IR-finiteness, interference terms between QED 
and QCD have to be taken into account, for both virtual and real (bremsstrahlung) 
contributions. Moreover, a new class of photon-induced tt production parton processes 
occurs, which for the LHC yield larger effects than the corrections to qq annihilation 
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and gg fusion. In size, the NLO QED contributions can reach the level of 5%. When 
combined with the rest of the EW corrections, the effects can become significantly large 
to find consideration for precision studies. 
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Figure 9: Differential cross sections (left) and relative correction 6 (right), as functions 
of the transverse momentum of the top quark (up) and of the parton energy (down), 
at the LHC, with no additional cuts. 
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Figure 10: Differential cross sections (left) and relative correction 5 (right), as functions 
of the transverse momentum of the top quark (up) and of the parton energy (down), 
at the Tevatron, with no additional cuts. 
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Figure 11: NLO QED contributions from the gg, qq and gj channels at the LHC for 
the pt and a/J distributions, including also cuts. 
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Figure 12: NLO QED contributions from the gg, qq and g'y channels at the Tevatron 
for the pt and \/I distributions, including also cuts. 
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Figure 13: Overall NLO QED effects in pp collisions at the LHC, for the distribution 
(up) and the a/5 distribution (down), without and with application of cuts. 
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Figure 14: Overall NLO QED effects in pp collisions at the Tevatron, for the p^ dis- 
tribution (up) and the a/J distribution (down), without and with application of cuts. 
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Appendix: Photon/gluon bremsstrahlung 



In the phase space shcing approach the phase space is divided into a region where 
the integrand is finite and into regions where the singularities occur. In the non- 
singular case the integration is performed numerically whereas in the singular regions 
the integration is carried out analytically in the approximation that the photon or 
photon-hke gluon is soft and/or coUinear to a charged fermion. 

We separate the soft and coUinear part of the singular regions by introducing two 
cutoff parameters AE and A^. In the soft part the photon/gluon energy satisfies the 
condition < AE <^ v^I, while in the coUinear part we have > AE and 9~^f < AO, 
where 6'-y/ is the angle between the photon and a charged fermion. 

In both regions the squared matrix elements for the radiative process factorize 
into the lowest-order matrix elements and universal factors containig the singularities. 
Thus, we can decompose the real corrections into 

dKeai = d^oft + da'^oii + dafi^ite, (22) 

where a = qq, gg. The coUinear contribution is zero in gg fusion, and in the case of qq 
annihilation only initial state radiation contributes since there are no mass singularities 
related to the final state bremsstrahlung. 

The soft part is combined with IR and mass singular virtual corrections to cancel the 
IR singularities proportional to log A and the mass singularities of double logarithms 
log^mg. The single logarithms logm^ are not compensated in the sum of virtual and 
real corrections and have to be handled by means of factorization. Single virtual and 
real contributions are dependent on the cutoff parameters AE and A^. However, the 
dependence has to cancel in the combination (22). 

The soft photon/gluon bremsstrahlung cross sections can be factorized into soft 
factors and the Born cross sections, 

da%, = ~daZrn-iC, + Q + C,,), 

d^foft - ~daZrn-Ct, (23) 

with Cq, Cf, Cqt referring to initial state radiation, final state radiation, and initial-final 
state radiation interference, respectively. 
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(24) 



with /3 = A^l — 4mf/s. The additional factor of 3 in the interference term comes from 
the gluon radiation contribution. 

The coUinear part of initial state radiation arises only for photons from the qq 
annihilation channel, and can be expressed as follows (see e.g. [45]), 



with the splitting function 
1 + 



In 



P,,{z) + (1 - ^) , (25) 



(26) 
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